Knight shifts around vacancies in the 2D Heisenberg model 
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The local response to a uniform field around vacancies in the two-dimensional (2D) spin-1/2 
Heisenberg antiferromagnet is determined by numerical quantum Monte Carlo simulations as a 
function of temperature. It is possible to separate the Knight shifts into uniform and staggered 
contributions on the lattice which are analyzed and understood in detail. The contributions show 
interesting long and short range behavior that may be of relevance in NMR and susceptibility 
measurements. For more than one impurity remarkable non-linear enhancement and cancellation 
effects take place. We predict that the Curie impurity susceptibility will be observable for a random 
impurity concentration even in the thermodynamic limit. 

PACS numbers: 75.10.Jm, 74.25. Nf, 75.20.Hr, 75.40.Mg 
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The deliberate use of substitutional impurities in 
strongly correlated electron systems has become a 
valuable tool for controlled studies of the underlying 
correlations 0, |M 1^ Q IS- Doping of non-magnetic Zn- 
ions in antiferromagnetic Cu02 planes is known to lead 
to a surprisingly large reduction of TcQ] and induces lo- 
cal staggered magnetic moments around the impurity 
sites 0- It has now become quite standard to analyze 
the local magnetic moments around static magnetic and 
non-magnetic impurities in an antiferromagnetic back- 
ground for a deeper understanding of the correlated 
states 11110. 

Theoretically the Knight shifts around impurities and 
vacancies have been studied for m any low-dimensional 
antiferromagnetsd S H H E [H [S , which typically 
show a strong enhancement of the antiferromagnetic or- 
der. The detailed behavior of the staggered magnetic 
moments can even give rather exotic results such as an in- 
creasing amplitude as a function of distanceQ and show 
renormalization effects [ill Il2j. The magnetization pat- 
tern can often be interpreted as the interference of in- 
coming and scattered quasiparticle excitations |llj or in 
terms of a pruned valence bond basisQ. 

The sum of all Knight shifts adds up to the total 
susceptibility xi- The resulting impurity susceptibility 
Ximp = Xi ~ Xq has been studied in detail in more re- 
cent theoretical studies HI [Illllllllllllllili^. In 
systems with long range magnetic order the leading tem- 
perature dependence of the impurity susceptibility from 
one single vacancy is given by a classical Curie spin which 
has been confirmed by numerical simulations for the 2D 
Heisenberg antiferromagnet [H Hsf. where a subleading 
log arithmic term has also been established [l^ 111 lisl lid 
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Here xo is the total susceptibility without any impurities, 
Ps is the spin stiffness, and the limit of large correlation 
length ^(T) and system size ^(T) > L ^ oo is assumed. 



We now want to examine how this impurity suscep- 
tibility is distributed on the lattice by considering the 
linear local response to a small uniform magnetic field 
(Knight shift) at each site 



(2) 



We find that long and short range patterns of the Knight 
shifts contribute to the impurity susceptibility differently. 
We are able to study the interference of several impuri- 
ties with non-trivial cancellation and enhancement effects 
over long distances, which allows us to make predictions 
for a finite impurity density. 

In general we can write as a sum of uniform and 
staggered parts on the lattice 



x(r)=Xu„i(r) + (-l)'^-+'-''Xstag(r), 



(3) 



the amplitudes of which are both slowly varying on 
the scale of one lattice spacing. In order to extract 
those two components we extrapolate the data on the 
even sublattice to the odd sublattice and vice versa 
and define Xuni/stag(r) = (Xovcn(r) ± Xodd(r))/2. The 
model is the 2D Heisenberg spin-1/2 antiferromagnet 
H ~ Jj2{ij)^i ' where (i, j) denotes nearest neigh- 
bor sites on a periodic square lattice. The quantum 
Monte Carlo program we developed uses the loop algo- 
rithm in a single cluster variety implemented in continu- 
ous time 21, 22, 23,,2^J, which gives efficient updates even 
at very low temperatures on a lattice of 100 x 100 sites. 
For convenience temperatures and energies are given in 
units of J = 1. 

The Knight shift data around one single vacancy is re- 
markably isotropic as a function of geometrical distance 
at all temperatures as is shown in Fig. ^ and the inset 
of Fig.|21for the staggered part. Slight anisotropic devia- 
tions near the edges (|r| 40 — 50) can be seen due to the 
finite system size and periodic boundary conditions. The 
remarkable isotropy on an anisotropic lattice is an indi- 
cation that the observed behavior is likely due to spin 
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FIG. 1: All data points of Xstag(r) as a function of distance |r| 
from one vacancy for T — 1, 0.5, 0.33, 0.2, 0.1 from below. In- 
set: All curves collapse to one almost universal, temperature 
independent shape after subtracting the limiting constant of 
the induced order Xoo for 0.02 < T < 0.2. 
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FIG. 2: (Color online) The uniform part 

Xuni(r) as a function of distance |r| for T — 
0.5,0.33,0.2,0.1,0.067,0.05,0.033,0.025,0.02 (from above at 
r = 10). Inset: The staggered part at T = 0.05 J 



wave excitations, which have an isotropic dispersion at 
low energies. At high temperatures T the impurity effect 
is confined within the correlation le ngth £(T) , which has 
an exponential dependence with 1/T 25,l2y,|23. Already 
for P = b the impurity affects the whole system in our 
case. This is commonly referred to as the zero temper- 
ature limit where Neel order effects can be observed in 
2D. In finite systems the effects of this order can be ob- 
served already for finite temperatures when the correla- 
tion length exceeds the system size ^ (T) > L (in our case 
when /? ^ 5). However, even at the lowest temperatures 
T we always remain in the thermodynamic limit, i.e. a 
large number of excited states contribute to a possible 
symmetry breaking and the energy level spacing remains 
very small compared to T. 

The origin of a classical Curie impurity susceptibility 
from a single vacancy in Eq. |^ can intuitively be un- 
derstood by considering the classical limit of an Ising 
system with N spins of size S (here S = ■^). Removing 
one spin on the even sublattice at (0,0) does not affect 
the long range order of such a system at all. In that 
case all the spins on the odd sublattice form one large 
effective spin of size NS/2, while the spins on the even 
sublattice have total spin of NS/2 — S. The difference of 
these two large spins will leave a net classical impurity 
spin of size S, which is free to rotate toward an applied 
magnetic field. In the 2D Heisenberg antiferromagnet a 
similar mechanism can be invoked, except that the long 
range order amounts only to about 61% of the classical 
limit HIH. As can be seen in Fig. we can indeed find 
a field-induced longitudinal long range staggered order 
which approaches a limiting constant Xoo = Xstag(oo)- 
The value of Xoo approaches a Curie behavior of about 
0.6 as can be seen in the inset of Fig. Q, consis- 
tent with the fact that only about 60% of the classically 



possible long range order couples to the field. The cor- 
responding contribution to the impurity susceptibility is 
given by 

Ximp-h- - 5](-l)^+^+^Xoo = Xoo 0.6^. (4) 

Since this contribution from long range order only 
amounts to 60% of the leading impurity susceptibility 
in Eq. the remaining 40% and subleading contribu- 
tions of Ximp must arise from a more local disturbance of 
Knight shifts in the vicinity around the vacancy. 

The uniform part of the Knight shifts is found to drop 
off very fast with the distance from the impurity to a 
value Xpurc which corresponds to the bulk susceptibility 
in a system without impurities xo/^^- This is shown 
in Fig. |5] for various temperatures as a function of dis- 
tance. Closer inspection shows that the behavior on dis- 
tance is not completely isotropic, since the enhancement 
is stronger along the lattice axes (at integer distances 
|r| = 1,2,3). The sum of this uniform part gives a net 
contribution to the impurity susceptibility 

Ximp-uni = ^ Xuni(r) " XO, (5) 

which turns out to contain the remaining 40% of the clas- 
sical Curie impurity susceptibility of l/12r in Eq. 
This is shown in Fig. |31 in comparison to other contri- 
butions to the impurity susceptibility, where Ximp-uni is 
seen to be proportional to /? and the sum Xoo + Ximp-uni 
follows the predicted form of the leading Curie contribu- 
tion 1/12T closely. 

Finally, there is a staggered pattern of Knight shifts 
localized around the impurity Xstag(r) — Xoo as shown 
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FIG. 3: The difFerent contributions to the impurity suscepti- 
bihty as a function of /3. The error in the uniform impurity 
susceptibility is about 0.2 because it is a difference of two 
large values. Inset: The ratio 12rxoo as a function of /3. 



in the inset of Fig. Interestingly, the dependence on 
temperature is very weak both in shape and magnitude. 
The net contribution to the impurity susceptibihty 



Ximp— stag 



(-1) 
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(Xstag(r) -Xoo) (6) 



is also small and contributes only to higher orders of Ximp 
as shown in Fig. |21 compared to the other contributions. 

Obviously Ximp = Ximp-stag + Ximp-lr + Ximp-uni- 

We now consider the case of two impurities that are 
separated sufficiently far away from each other, so that 
any quantum interference can be neglected. As long as 
the correlation length ^ is smaller than the inter-impurity 
distance no interference can be found and the impurity 
contributions to the Knight shifts and the susceptibility 
simply add. In the low temperature limit the correlation 
length grows beyond the system size and the staggered 
order extends through the entire system. In this case it 
makes a significant difference if the two vacancies are on 
the same or on different sublattices. 

For two vacancies on different sublattices the field- 
induced long range order cancels exactly, giving no con- 
tribution to the impurity susceptibility. Remarkably, the 
uniform contribution is also strongly affected as shown 
in Fig. 2] for impurities located at (0, 0) and (48, 49) in 
comparison to the uniform Knight shift data of just one 
impurity. Even though the impurities are located a dis- 
tance of |r| ~ 69 sites apart, we find a strong reduction 
of the uniform Knight shifts within just a few lattice sites 
around each impurity in Fig. 0] This reduction goes far 
beyond simple interference and must be taken as a col- 
lective effect of the entire system. The resulting impurity 
contribution Ximp-uni from the uniform Knight shifts in 
Eq. ISJ is small and has no Curie-like temperature be- 
havior, i.e. only the subleading corrections remain. On 
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FIG. 4: (Color online) The uniform part of the Knight shifts 
around an impurity at (0, 0) as a function of distance at 
(3 = 20. For a second impurity on the same sublattice (48, 48) 
(dotted) we find twice the uniform enhancement compared to 
the single impurity case (dashed). For a second impurity on 
the other sublattice (48, 49) (solid) there is a strong reduc- 
tion. Inset: Staggered part for two impurities on different 
sublattices. 



the other hand, the local staggered enhancement of the 
Knight shifts is remarkably unaffected by interference ef- 
fects, but has the opposite phase around each impurity as 
shown in the inset of Fig. 21 Therefore, the sum in Eq. ® 
gives exactly twice the magnitude of Ximp-stag compared 
to one single impurity (the opposite phase cancels with 
the minus sign from the different sublattice in Eq. ©). 
This part Ximp-stag now becomes the main contribution 
to the impurity susceptibility, but has no divergent Curie- 
like behavior as discussed above. This is in agreement 
with findings in Ref. that the impurity susceptibility 
of two impurities on different sublattices corresponds to 
subleading terms only. 

For two vacancies on the same sublattice we find a 
doubling of both the induced longitudinal long range 
order 2xoo and the uniform enhancement as shown in 
Fig. ^ Since two vacancy sites contribute in Eqs. (0)) 
and (jnj this implies that the Curie-like impurity sus- 
ceptibilities Ximp-lr and Ximp-uni are quadrupled as if 
two virtual spins form a triplet. On the other hand, 
the local staggered enhancement of Knight shifts close 
to each impurity again remains almost completely unaf- 
fected, compared to the single impurity case in the in- 
set of Fig. n The Curie-like contributions Ximp-ir and 
Ximp-uni therefore increase quadratically with the num- 
ber of same-sublattice impurities N, while the subleading 
contribution Ximp-stag increases linearly with N indepen- 
dent of sublattice Ximp N^/12T + NO (ln(C/T)). 

In the case of two impurities that are very close to each 
other, quantum effects are bound to also play a role and 
a comparison with the interference of two single impu- 
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rities is not useful. For two neighboring impurities we 
find that also subleadin g co ntributions to the impurity 
susceptibility are smaller [isj. 

In an infinite system with a small random impurity 
concentration p the impurities act independently at high 
and intermediate temperatures (assuming an average dis- 
tance of at least several lattice sites). As the tempera- 
ture is lowered, however, vacancies within the distance 
of the exponentially growing correlation length ^ start 
to feel each other, i.e. the effective number of corre- 
lated impurities within an antiferromagnetic domain is 
N cx ^^p. The subleading contribution from the lo- 
cal staggered part Ximp-stag always adds independent 
of sublattice and will therefore be of order N, but is 
not divergent with temperature. For the Curie-like con- 
tributions Ximp-ir and Ximp-uni the difference of impu- 
rities on the two sublattices plays an important role, 
which is of order \Na — Nb\ ^ V^Q- Our simula- 
tions show that the induced long-range order is pro- 
portional to the excess of impurities on one sublattice 
\Na — Nb\Xoc £,s/pXoo while the uniform enhancement 
is proportional to ±{Na — Nb) on the sublattice A/B up 
to subleading terms, i.e. negative on the minority sublat- 
tice. The contributions Ximp-ir and Ximp-uni in the sums 
Q and lO therefore both increase quadratically with 
the difference {Na — NbY ^ Hence we predict that 
a Curie susceptibility of pS^ /3kBT survives even in the 
thermodynamic limit, contrary to conjectures made in 
Ref. In NMR experiments induced staggered Knight 
shifts of order ^^Xoa^ throughout the lattice should be 
directly observable. 

In conclusion we have analyzed the local response 
around vacancies in the 2D Heisenberg antiferromagnet 
in the low temperature limit. These Knight shifts can 
be analyzed in terms of staggered and uniform contri- 
butions, which are mostly isotropic on the 2D lattice. In 
the presence of a small magnetic field a single vacancy in- 
duces a staggered pattern of Knight shifts over the entire 
lattice in the low temperature limit, giving a longitudinal 
Neel order ~ 0.6/12T. This field-induced order cor- 
responds to a free classical moment accounting for 60% of 
the Curie impurity susceptibility. The remaining 40% of 
the magnetic moment is found to be located in a uniform 
enhancement of Knight shifts in the immediate neighbor- 
hood of the vacancy. Additionally a staggered pattern 
in a finite range around the vacancy is also found, which 
turns out to have a largely temperature independent, uni- 
versal shape and gives rise to subleading corrections of 

Ximp- 

When more than one vacancy is introduced in the low 
temperature limit, the Curie terms in the impurity sus- 
ceptibility cancel if the vacancies are on different sub- 
lattices with a remarkable long-distance reduction of the 
uniform Knight shifts in the close vicinity of each va- 
cancy. For vacancies on the same sublattice both the 
induced long range order and the uniform part increase. 



leading to a Curie-like impurity susceptibility that in- 
creases quadratically with the number of impurities. Both 
uniform and long range staggered parts appear to be con- 
tributing to virtual classical impurity spins that are ferro- 
magnetically/antiferromagnetically linked for impurities 
on the same/opposite sublattice. The local staggered 
enhancement around each impurity remains unaffected 
by interference and temperature in any case, contribut- 
ing to subleading corrections of Ximp- In the thermody- 
namic limit cancellation and non-linear enhancement ef- 
fects even out from statistical arguments, leaving a Curie 
susceptibility of p/12T. 

Some of the results can also be generalized for other 
quantum antiferromagnets, like Heisenberg models in 
higher dimensions and/or with higher spin. In partic- 
ular, a vacancy will always induce a longitudinal anti- 
ferromagnetic response Xoo over the entire lattice, the 
magnitude of which is linked to the order parameter 
as T ^ (80% in Eq. (gj) for a 3D spin-1/2 model). 
The remaining impurity susceptibility in Eq. must be 
found in a uniform enhancement around the vacancies. 
The virtual classical impurity spins are antiferromagnet- 
ically/ferromagnetically linked depending on the sublat- 
tice as seen above. In the spin-1/2 chain no long-range 
order and no Curie impurity susceptibility is found. How- 
ever, the induced Knight shifts increase with a distance 
from the vacancy and boundary terms give a Curie 
contribution that is logarithmically reduced for low tem- 
peratures from non-universal boundary effects [2^. For 
vacancies in quantum antiferromagnets with a gap, a 
more local antiferromagnetic pattern is induced [a, 
leading to a Curie susceptibility corresponding to vir- 
tual quantum spins, which arc not linked to each other. 
The effect of substitutional impurities with higher spin 
is much less clear, although a coupling to the antiferro- 
magnetic order can again be expected. 
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